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Abstract 

We report results of diffusion Monte Carlo calculations for both '^He ab- 
sorbed in a narrow single walled carbon nanotube {R = 3.42 A) and strictly 
one dimensional ^He. Inside the tube, the binding energy of liquid ^He is 
approximately three times larger than on planar graphite. At low linear den- 
sities, ^He in a nanotube is an experimental realization of a one-dimensional 
quantum fluid. However, when the density increases the structural and en- 
ergetic properties of both systems differ. At high density, a quasi-continuous 

liquid-solid phase transition is observed in both cases. 
PACS numbers:05.30.Jp,67.40.Kh 
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Since their discovery by Ijima [||] in 1991, carbon nanotubes have received a great deal 
of attention. Basically, they are the result of the seamless rolling up of one or several 
graphite sheets over themselves Depending on the relative orientation of the rolling 

axis with respect underlying graphite structure, one can have different types of nanotubes 
0: armchair, zig-zag and chiral with different radii and different mechanical and electrical 
properties. Nowadays, it is possible to obtain high yields of nanotubes (single and multiple 
walled), with a variety of diameters ranging from 7 to 40 A [|] and lengths up to ~ 1000 
times larger. 

One of the most attractive features of carbon nanotubes is the possibility of filling with 
different materials both their inner cavities and the interstitial channels among them 
The interest in this field is twofold. On one hand, the expected increase in the particle- 
substrate potential energy with respect to a fiat carbon surface has suggested the use of 
nanotubes as storage devices for molecular hydrogen in fuel cells [§,0]- On the other, more 
theoretical, nanotubes provide a reliable realization of one- dimensional systems in the same 
way that a substance adsorbed on graphite manifests trends that are characteristic of a two- 
dimensional medium. If the nanotubes are filled with light atoms (He) or molecules (H2) and 
the temperature is low enough, one is dealing with quasi-one dimensional quantum fluids. 
Such an experimental realization has been carried out for the first time by Yano et al. |TlI 



in a honeycomb of FSM-16. This is a mesoporous substrate with tubes approximately 18 A 
in diameter. Using a torsional oscillator, this group proved the existence of superfluidity of 
the ^He atoms absorbed in the pores below a critical temperature of ~ 0.7 K. More recently, 
Teizer et al [0 have studied experimentally the desorption of ^He previously absorbed in the 



interstitial sites of carbon nanotube bundles. In this case, the data points unambiguously 
to the one-dimensional nature of the helium inside the nanotubes. 

From a theoretical point of view, it has been recently established using both the hyper- 
netted chain (HNC) variational approach [|I^ and the DMC method ||14[ that strictly one 
dimensional (ID) ^He is a self-bound liquid at zero temperature. However, contrary to the 
situation for dilute classical gases [|T^ -|l6| , there are no many body calculations of quantum 



fluids inside nanotubes yet. In this work, we address the question of the quasi-one dimen- 
sionahty of ^He absorbed in a tube by a direct comparison between the results of diffusion 
Monte Carlo (DMC) for strictly ID ^He and ^He inside a nanotube of radius equal to 3.42 
A, which corresponds to a (5,5) armchair tube in the standard nomenclature 0. 

The DMC method [0,0 solves stochastically the A^-body Schrodinger equation giving 
results that are exact for bosonic systems as liquid ^He, provided that the interatomic po- 
tential is known. In the present calculation, we have used the HFD-B(HE) Aziz potential 
for the He-He pair interaction [^, and the potential given by Stan and Cole in their 



study of Lennard- Jones fluids in tubes for the He-tube one. Basically, they consider the 
nanotubes as smooth cylinders by making a 2;-average of the corresponding sum of all the 
C-He interactions. Thus, the potential felt by a particle only depends on its distance to the 
center of the cylinder. This is a simplification, but one would expect the error involved to be 
small since the helium atoms are much larger than the C-C distance. In fact, the differences 
in energy and position between a ^He atom in the smooth cylinder model and the same 
particle considering its interaction with the surrounding individual carbons are about 1% 



for the tube considered here IE 



The efficiency of the DMC method is greatly enhanced by introducing a trial wave 
function ^^(R) that acts as an importance sampling auxiliary function. In ID '^He we have 
used a two-body Jastrow wave function 



vpi°(R) = v[/j(R) (1) 



with ^'j(R) = rii^ exp 



2 \njj 



, whereas liquid ^He inside nanotubes requires the addi- 



tional introduction of a one-body term 

*T(R) = ^j(R)^e(R) (2) 

with ^'c(R) = nf^ exp(— erf ) (r, being the radial distance of the particle to the center), that 
accounts for the hard core of the helium-nanotube interaction. Using the variational Monte 
Carlo method (VMC) we have optimized the parameters b and c at low densities around 



the equilibrium. The values obtained, b = 3.067 A and c= 2.679 A~^, show a negligible 
dependence with the density and therefore they have been used everywhere in the DMC 
calculations. In all the simulations = 30 atoms have been used, a number that has 
proved to be large enough to reduce the size effects to the level of the statistical errors 
reported. 

The variational HNC equation of Krotscheck and Miller |13| on ID liquid ^He points to 
the existence of a liquid-solid phase transition at high density. We have explored this feature, 
that is only possible at zero temperature, by using a solid trial function which results from the 
product of ^^°(R) and ^"^(R) by a z-localized factor ^s(-R) = exp{-a{zi - Zisf). The 
solid sites Zis are equally-spaced points along the z direction which is both the longitudinal 
axes of the tube and the line of the ID system. In both solid systems, a VMC optimization 
at high densities leads to values b= 2.939 A and a = 0.612 A~^, with c= 2.908 A~^ in the 
tube case. 

The energy per helium atom versus the linear coverage. A, for the ID (open squares, 
energy scale on the right) and the tube (full squares, energy scale on the left) is shown 
in Fig. 1. The two curves have have been drawn for the full square with the lowest A to 
coincide with the open symbol for the same He density. One observes that for A < 0.05 
A-^ both curves are similar, but for larger concentrations the tube curve is located below 
the other one. A similar phenomenology appears in the comparison between the energies of 
purely 2D ^He and ^He adsorbed in graphite. As for this system, the difference in energy 
between ^He in a nanotube and ID ^He is always negative with an absolute value that 
increases with the density. Both in graphite and in nanotubes this increase with respect to 
the 2D and ID systems is mainly due to the emergence of their actual 3D nature. Beyond 
this qualitative agreement between ^He adsorbed on graphite and inside carbon nanotubes, 
there are significant differences in the values of the binding energies in the two systems. 
The binding energy of a single "^He atom in graphite is E^= 140.74 K |^, whereas in 
the nanotube we are studying is roughly three times larger, E]^ = 429.97 K, a significant 
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difference that has been dramatically observed in the desorption experiment of Teizel et 
al |T^. On the other hand, the departure of the real 3D systems (^He in graphite or a 
nanotube) from the idealized 2D or ID liquids can be quantified by means of the parameter 

^T(G) ^ {E'^^^'^ -EI'-^^) -^i(2)D 
(^T(G) _^T(G)^ 

where T stands for the tube and G for the graphite adsorbents and E is the energy per 
particle in the system under consideration. Around the respective equilibrium densities one 
obtains A"^ = 90 % and = 6 %, a large difference that indicates that the ID representation 
of ^He inside the nanotube is worse than the 2D modelization of ^He in planar graphite. 

Up to A = 0.15 the energies per particle (e = E/N) of both ID ^He and ^He inside 
the tube may be well fitted by a third-degree polynomial 

The optimal values for the parameters A, B, Aq, and Cq are reported in Table I. The linear 
equilibrium densities Aq of both systems are close, Aq° = 0.062 and Aq = 0.079 
whereas the energy difference between A = and A = Aq is significantly different, 0.00364 
and 0.018 K for ID ^He and ^He inside the tube, respectively. The latter results point 
again to a large enhancement of the binding energy of ^He inside the tube with respect to 
the ID system. On the other hand, the equilibrium density of liquid ^He inside the tube 
(po = 0.0022 A^^) is much smaller than the one in homogeneous 3D liquid ^He (po = 0.022 

A-3). 



In agreement with the DMC calculation of Stan et al. |14] and the variational one of 



Krotscheck and Miller |T3|, ^He selfbounds in a ID array but with a binding energy 



0.0036±0.0002 K) much smaller than that in 2D (-0.897±0.002 K) @ and 3D (-7.267±0.013 
K) [0. It is worth noting that such a small total energy results from a big cancellation 
between the potential and kinetic energies. At Aq, we have T/N = 0.2706±0.0004 K and 
V/N = -0.2742±0.0004 K. In fact, the influence of the ^He interatomic potential in this 
system is very large. A calculation at the equilibrium density Aq for the ID system using 



the HFDHE2 Aziz potential p3[ indicates that ^He is still a liquid, but the total energy is 



a factor two smaller (-0.0018 ± 0.0003 K, with a potential energy -0.2724±0.0004 K and 
the same kinetic energy). This sizeable differences partially explain the discrepancies of the 



DMC calculation of Stan et al. and ours with the results of Krotscheck and Miller |T3 
who used the HFDHE2 Aziz potential. 

From the values of the energy, one can obtain the linear system pressure, p\ = X'^de/dX, 
and estimate the same property for helium inside the cylinder as p = px/nR^. Fig. 2 
displays this observable as a function of the ^He density. The range corresponds to a liquid 
structure both in the ID and tube cases (see discussion below). One can see that the 
pressure increases faster in a pure linear arrangement of atoms than in a tube. This can be 
understood if one considers that in a narrow nanotube it is possible to avoid the repulsive 
core of the nearest neighbors by shifting transversely the helium positions, a situation that is 
obviously not possible in ID. Also interesting is the comparison between the sound velocity, 
c(A) = [l/m((9P/9A)]^/^ at their respective equilibrium densities. The values are Cid = 7.98 
± 0.07 m/s and ct = 14.2 ± 0.8 m/s, in both tiny fraction of the corresponding 

2D (c = 92.8 m/s) ^ and 3D (238.3 m/s) |g ^He hquids. The spinodal point can be 



obtained as the density at which the speed of sound becomes zero. According to our results 
the spinodal points are located at Am = 0.047 ± 0.001 A^^ and At = 0.059 ± 0.001 A'^ 

Another aspect that has deserved our attention has been the existence of a liquid-solid 
phase transition at high densities. Evidences of this phase transition, that is only possible at 
zero temperature, appear in a variational calculation of ID ^He [|13|. A comparison between 
the DMC energies for the liquid and solid phases is given in Table II. One can see that in 
both systems, the energy per particle when localization is imposed (a ^ 0) is below the 
corresponding to a liquid structure (a = 0) for lineal densities greater than 0.358 A^^. By 
means of the Maxwell double tangent construction, one would be able in principle to tell the 
solid from the liquid and to obtain the freezing and melting densities. Unfortunately, the 
energy differences between the 2;-localized and the liquid structures are too small to allow us 
to carry out a meaningful calculation. Our results indicate that for large enough densities 
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(A > 0.358 A-i) both the ID and the tube arrangements have a locahzed (sohd) phase 
and that the discontinuity in the density (if any) is surely very small. This is in agreement 
with the results discussed by Withlock et al. ||2lj| about the reduction of the size of this 
discontinuity from three to two dimensions: in 3D is fairly large, being considerably smaller 
in a purely 2D system. In ID, we observe a further reduction towards a continuous or a 
quasi-continuous transition. It is also remarkable that ^He inside the carbon tube remains 
a liquid up to a much larger pressure (around 5 times) than in bulk liquid ''He(~ 2.6 MPa). 

Information on the spatial distribution of the ^He atoms may be drawn from the two-body 
radial distribution function along the z direction, gz{r). The functions gz{r) for ID ^He and 
^He in the tube are shown in Fig. 3 at several linear densities. Near the equilibrium density 
(A = 0.08 A^^, lower part of the figure) gl^{r) is quite similar to gj{r), as corresponds to 
a quasi-one dimensional system. The same could be said in a broad range of densities, as 
it can be seen in the curves for A = 0.182 A^^ (middle part of the figure). On the other 
hand, in the solid phase (A = 0.406 A"^), gl^{r) and gj{r) are different: in this case, the 
3D nature of ^He inside the tubes produces a significant decrease in the localization with 
respect to the ID result. 

In conclusion, we have compared the properties of strictly ID ^He with ^He inside a nar- 
row carbon nanotube using the diffusion Monte Carlo method. For a wide range of densities, 
''He is a liquid in both systems, and also in both cases a quasi-continuous liquid-solid phase 
transition has been observed. In accordance with recent experimental determinations, the 
present calculation evidences a quasi-one dimensional behaviour of '^He inside a nanotube 
but significant differences with the ideal ID system appear, specially when the linear den- 
sity is increased. The origin of these differences is mainly the existence of the additional 
transverse degree of freedom that helium atoms have inside a nanotube. 

One of us (M.C.G.) thanks the Spanish Ministry of Education and Culture (MEC) for 
a postgraduate contract. This work has been partially supported by DGES (Spain) Grant 
No. PB96-0170-C03-02. 
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TABLES 



Parameter ID ^He ^He in a tube 

Ao (A-i) 0.062 ± 0.001 0.079 ± 0.003 

eo (K) -0.0036 ± 0.0002 -429.984 ± 0.001 

A (K) 0.0156 ± 0.0009 0.048 ± 0.006 

B (K) 0.0121 ± 0.0008 0.0296 ± 0.009 

x'^/f 2.2 0.24 
TABLE I. Parameters of Eq. 4 for the two systems studied. 



A (A-i) 


E/N (ID, a = 0) 


E/A (ID, a / 0) 


£;/A (T, = 0) 


E/N (T, o 7^ 0) 


0.406 


123.726 ± 0.012 


123.561 ± 0.012 


-350.155 ± 0.030 


-350.20 ± 0.02 


0.380 


67.070 ± 0.011 


67.000 ± 0.009 


-382.282 ± 0.016 


-382.321 ± 0.012 


0.358 


37.602 lb 0.008 


37.596 ± 0.007 


-401.873 ± 0.013 


-401.844 ± 0.010 


0.338 


21.881 ± 0.007 


21.904 ± 0.005 


-413.091 ± 0.014 


-413.061 ± 0.012 


0.320 


13.240 ± 0.005 


13.258 ± 0.006 


-419.551 ± 0.011 


-419.493 ± 0.010 



TABLE IL Energies per particle at large A for the systems studied. All the energies are in K. 
See text for further details 
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Figure captions 

1. Energy per particle {E/N) versus the linear concentration (A), for the two systems 
we have studied: a strictly one dimensional system (open squares, right energy scale), and 
a (5,5) armchair tube (full squares, left energy scale). In the first case, the error bars are 
less than the size of the symbols. Both energy scales are in K. 

2. Pressure at high helium densities for both systems (one dimensional, dashed line, left 
scale; nanotube, full line, right scale). 

3. Pair distribution function along the z coordinate, g'^^^^\r). Full lines correspond to 
the narrow tube at 0.08 A'^ (bottom), 0.182 A'^ (middle) and 0.406 A"^ (top), and dashed 
lines to the purely linear system at the same densities. 
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